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ABSTRACT 

We have studied the photometric properties of four fields around the high- 
redshift quasar pairs QP1310+0007, QP1355-0032, QP0110-0219, and QP0114- 
3140 at z ~ 1 with the aim of identifying large-scale structures- galaxy clusters 
or groups- around them. This sample was observed with GMOS in Gemini North 
and South telescopes in the g', r', i', and z' bands, and our photometry is complete 
to a limiting magnitude of i' ~ 24 mag (corresponding to ~ M*, + 2 at the 
redshift of the pairs). Our analysis reveals that QP0110-0219 shows very strong 
and QP1310+0007 and QP1355-0032 show some evidence for the presence of 
rich galaxy clusters in direct vicinity of the pairs. On the other hand, QP0114- 
3140 could be an isolated pair in a poor environment. This work suggest that 
z ~ 1 quasar pairs are excellent tracers of high density environments and this 
same technique may be useful to find clusters at higher redshifts. 

Subject headings: galaxies: clusters: general — galaxies: distances and redshifts 
- galaxies: high-redshift — large-scale structure of universe — quasars: general 
— X-ray: galaxies: clusters 



Introduction 



The study of galaxy populations in high-redshift large-scale structures can give us im- 
portant clu^s_jboui_the star in such environ- 
ment s J White fc Frenklll99ll : lBekkilll998l : kodama et allll998l : Ivan Dokkuml 1200,4 iMei et al 
20061 ). However, the detection of distant galaxy clusters is not trivial. 



1 Departamento de Astronomia, Instituto de Astronomia, Geoffsica e Ciencias Atmosfericas, Univer- 
sidade de Sao Paulo, Rua do Matao 1226, Cidade Universitaria 05508-900, Sao Paulo, SP, Brazil; na- 
talia@astro.iag.usp.br 

2 Department of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, 
United Kingdom 

3 Gemini Observatory/ AURA, Southern Operations Center, Casilla 603, La Serena, Chile 



- 2 - 



At high redshifts, the use of techniques like the red-cluster sequence (IGladders fc Yee 



20051 ). is less efficient. The reason is the proportional decrease of the number of red galaxies 
in clusters for increasing redshift (e.g. iButcher fc Oemlerl Il984l ). Therefore, in order to 
detect high-z clusters through optical imaging, one requires at least photometric redshift 
information, which implies imaging in four or more bands. Additionally, in order to avoid 
the high observational cost of observing large areas of the sky, one can use several indicators 
of the presence of high-z clusters to select the fields to be observed. 



Several of these tracers, like extended X-ray emission (e.g., iRomer et al.l l200ll ). the 
Sunyaev-Zeldovich decrement or bright radio-emitting galaxies have been largely used to 
trace clustering of galaxies. These techniques have strengths and drawbacks. The first two, 
for example, depend on the presence of a hot intra-cluster medium, which may bias samples 
against recently forming clusters. In this paper, we will consider another possible tracer of 
the presence of clusters: physically close pairs of quasars. 

Quasars are relatively rare astronomical objects and hence, if they are distributed fol- 
lowing galaxies, the presence of two or more such objects in a relatively small volume should 
be a good indicator of a rich environment. Actually, in structure form ation scenarios with 
bias between barionic and dark matter distribution (je.g.. Kaiserl 119841 ) it is expected that 
high redshift objects form in large high-redshift density fluctuations and, therefore, such 
correlation between quasar concentration and clusters is somewhat expected, unless, for 
some reason, quasars avoid clusters. However, most observational evidence shows that high 
redshift quasars do tend to follow the overall large scale structures. 

Whet her q uasar s inhab it or not high density regions in low redshifts is a subject of dis- 
pute. IColdwell et al.l (120021 ). for example, claim that at 0.1 < z < 0.25, quasars (both radio 



l oud a nd radio-quiet) tend to reside in low density regions. On the other hand lMullis et al. 



( 120041 ). using a sample of X-ray s elected quasars, concl ude that those objects trace closely 
the underlying mass distribution. ISochting et al.l (120021 ) also points out that 0.2 < z < 0.3 
quasars follow the large-scale structure traced by galaxy clusters, but they also note the 
complete absence of radio-quiet QSO's at the very center of galaxy clusters. 

At higher reds hift, however, mos t observational re s ults s uggest that quasars prefer 
groups or clusters (IHall fc Greenl 19981: IWold et al.l |2000| . 120011 ). One very convincing ex- 



Haines et a 



ample is the structure found by 
quasar belonging to a large quasar structure ( 



(200 1J) at z = 1.226 around a radio-quiet 



Clowes fc Campusanol Il99ll . |l994j) . The same 



behaviour appears to be followed by radio-loud quasars. A good example is the work by 
Sanchez fc Gonzalez- Serrano! (120021 ) , who found a h i ghly s ignificant excess of galaxies around 
radio- loud quasars at 1.0 < z < 1.6. iTanaka et al.l (120011 ) also points in the same direction 
by reporting an overdensity of galaxies around a quasar concentration at z ~ 1.1. An excep- 
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tion is the work by ICoil et al.l (120071 ) who, through an analysis of the clustering of quasars 
and galaxies at 0.7 < z < 1.4, concluded that quasars and blue galaxies are found in the 
same environment, whic h differs from that o c cupie d by the red galaxy population. Regarding 
specifically quasar pairs, IZhdanov fc Surdejl (120011 ) found an statistically significant excess of 
high-redshift quasar pairs with separations between 1 and 5 Mpc in projected distance. This 
suggests that such quasar pairs belong to sizable physical structures (precursors of today's 
clusters and superclusters of galaxies) and therefore, they c an be used as tracers o f high- 
redshift large-scale structures. Going to even larger redshifts, iDjorgovski et al.l (120031 ) found 
that a quasar pair at z = 4.96 is associated with a large-scale structure. Thus, an interesting 
form to search for high-redshift clusters and other large-scale structures is examining the 
environment inhabited by quasar pairs. 

In this work, we describe a multi-color photometric study of the field around four quasar 
pairs at z ~ 1, using the instrument GMOS in both Gemini North and South telescopes. 



One of the pairs in our sample, QP0110-0219, has been previously studied by lSurdej et al 



(119861 ). who found hints of the presence of a cluster around it. The new data we present 
here allow us to confirm this claim. There are no studies in the literature for the other three 
quasar pairs. 

The outline of this paper is the following: in Section 2 we describe the sample and 
the data reduction procedures. The galaxy photometry is discussed in Section 3. Section 4 
outlines our approach to obtain photometric redshifts and presents then application to our 
galaxy sample. The environments of the quasar pairs are discussed in section 5. Finally, in 
Section 6 we summarize our results. Throughout this paper, we adopt a ACDM concordance 
cosmology with Q m = 0.3 and = 0.7, and we use the value h = 0.7 in the Hubble constant, 
H = 100 h km s" 1 Mpc" 1 . 



2. Observations and data reduction 
2.1. Sample selection 

In this pa per, we study a sample o f four fields around quasar pairs at z ~ 1. We selected 



the pairs from IVeron-Cetty fc Veronl (120011 ) quasar catalog considering redshift differences 
smaller than 0.01 and projected angular separations smaller than 300 arcsec. We did not 
consider pairs with angular separations smaller than 15 arcsecs to avoid including gravita- 
tional lens. With these parameters, we found 84 quasar pairs. Five of them had redshifts 
between 0.9 and 1.0, and four were observed with Gemini telescopes. The main sample char- 
acteristics are shown in Table [TJ It includes: the quasar names, their coordinates, redshifts, 
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angular separation, and the name adopted for the pairs in this paper. 

We have checked the spectra of the quasars in our sample to certify that we indeed 
did not pick any cases of gravit ationally len s ed im ages of one only quasar. The parity 



of QP0110-0219 is discussed by ISurdej et al.l (119861 ) . Considering the redshift difference 



and spectral characteristics, they conclude that the quasars Q 0107-0235 and PB 6291 are 
different objects. For the o ther pairs, quasar spectra are available in the 2dF QSO Redshift 
Survey ( Croom et al. 2004 1^1. Our visual examination of the spectra indicates that also in 



this case the differences in redshifts and spectral characteristics suggest that they are indeed 
different objects and not lensed images of the same quasar. Moreover, in QP1310+0007, 
QP1355-0032, and QP0110-0219 one of the quasars is radio-loud and the other is radio- 
quiet. Consequently, we are confident that none of the pairs in our sample are produced by 
gravitational lensing. 



2.2. Imaging and data reduction 

The four fields in Table [1] were observed with GMOS N and S mounted on Gemin i 



telescopes. The imaging was done in four filters of the SDSS system ( IFukugita et al.lll996l ): 
g', r', i', and z' . The log of observations is presented in Tabled which shows the telescope 
used, the exposure time, and the Gemini program identification number. All observations 
were performed in photometric conditions. The typical FWHM for point sources was ~ 0.7 
arcsec in all images. 

Data reduction was performed using the Gemini IRAF __ package. The images were 
bias corrected, flat fielded, and fringe corrected in the standard way. After that, they were 
combined and cleaned of cosmic ray events and bad pixels producing, then, the final images, 
appropriate for science analysis. 



3. Photometry and object detection 

We have used the IRAF package daophot to calculate the photometric zero-point for 
each band in each field in the AB SDSS photometric system. The calibration was made 



x http: / /www. 2dfquasar.org/Spec_Cat/2qzsearch2. html 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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using stars from the Landolt catalog (lLandortlll992l ) also calibrated in the AB SDSS system. 
Using the dispersion in the magnitudes of the stars, we have estimated the accuracy of the 
magnitude zero-point as 0.01 in g' and r' bands, 0.02 in i' band, and 0.03 in z' band. 



We have used SExtractor (IBertm fc Arnautslll996l ) to detect objects over the final image 
frames. First, we ran the program on the images of each photometric band and selected the 
image that showed the highest number of detected objects. Second, using such image as 
reference, we ran the program again in "dual image mode". We used a top-hat filter and 
detected objects above 1.5 a, which corresponded to median isophotal levels of 27.1, 26.4, 
26.4, and 25.4 mag arcsec -2 in g', r', i', and z', respectively. In order to run the program 
in " dual image mode" , it was necessary to align the images. Thus, because of rotations and 
shifts, parts of the images near the borders were lost. 



Positions and magnitudes (total and aperture) were obtained for all objects present 



in all bands for each pair. We adopted 3 arcsec aperture magnitudes, m ap , to measure 
colors. Aperture magnitudes were also obtined to compare our data with others in the 
literature. For the total magnitude of an object, we have addopted a color corrected isophotal 
magnitude. For example, if the objects were detected in the g' band, then g' = g iso , r' = 
9iso — {gap — r a p), etc., where rrii SO is the isophotal magnitude given by SExtractor. After 
measuring the magnitudes th ey were corrected for Galaxy extinction, with the absorption 
coefficients A\ obtained from ISchlegel et al.l (119981 ) using NED and interpolated to GMOS 
bands. 

We have used the class-star parameter of SExtractor, which ranges from (galaxies) 
to 1 (stars), to separate stars from galaxies. Figure [T] shows this parameter versus the i' 
magnitude for the pair QP0110-0219. The star symbol in the plot represents all objects with 
FWHM < seeing. If we consider a ll objects with class-star < 0.8 as galax ies, a threshold 
often adopted in the literature (e.g. iKodama et al.l 12004 ; ICaputi et al.ll2006l ). we find that 2 
% of the objects with FWHM < seeing have class-star < 0.8 and that 7 % of the objects 
with FWHM > seeing have class-star > 0.8. On the other hand, if we adopt a threshold 
of 0.9, we have a similar contamination (~2 %) of the galaxy sample and only 4 % of the 
objects with FWHM > seeing have now class-star > 0.9. The results for the other fields are 
similar. We then decided to adopt the class-star value of 0.9 to separate stars from galaxies. 
The same criteria was adopted by ICapak et al.l (120041 ) to determine number counts in the 
HHDFN. 

In order to estimate the completeness magnitude of the observations, we have ploted 
the logarithmic number of detected objects as a function of the total magnitude in the band 
used for detection. From visual inspection of the turnover magnitudes, we estimated that 
the observations are complete down to %' = 24 for QP1310+0007 and QP0110-0219, g' = 25 
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for QP1355-0032, and g' = 24.5 for QP0114-3140. 



It is interesting to know how the magnitudes above compare with those of a M*-galaxy 
at the redshi f t of t he pairs. We have estimated the value M* in two ways, as follows. 
Ellis fc Joned (120041 ) obtained K* ~ 18 for clusters of galaxi es with redshifts betw een 0.8 
and 1.0. Considering the value for (I — K) ~ 2.9 obtained by IStanford et al.l (120021 ) for the 
cluster 3C 184 (z = 0.996), we hav e I* ~ 20.9. On the other hand, using spectrophotometric 
synthesis models, iFukugita et al.l (119951 ) obtained (i' — I c ) ~ 0.7 for galaxies at z = 0.8. 
Then, we obtain i'* ~ 21.6. I n the second case, we may consider the Coma cluster as 
representative of a z = cluster. iMobasher et al.l (120031 ) st udied its lumino si ty fun ction and 
found Mft ~ —21.79 + 5 log h e5 . For galaxies at z — 0, IFukugita et al.l (119951 ) obtained 
(r' — R c ) oj 0.22 and (r' — i') ~ 0.30, therefore M*, ~ —21.71 for the cosmology adopted here. 
We have calculated i'* with 



m 



M(z = 0) + 5 log d L [Mpc] + 25 + k(z) + e(z) 



X 



where dt is the luminosity distance, k(z) is the k-correction and e(z) is the evolution cor- 
rection. In the cos mology adopted here, ~ 6000MjP c at z ~ 1. Using k(z) and e(z) 
values published by IFukugita et al.l (119951 ) and iPoggiantil (119971 ) respectively {kit ~ 0.9 and 
tii ~ —1.3) and the value obtained to M*, at z = 0, we obtain i'* ~ 21.8. Considering the 
uncertainties of the approaches, the agreement of the two values is very goo d. We have then 
adopt ed the mean value %'* ~ 21.7. This value is similar to that obtained by lBlakeslee et al. 



( 120061 ) (ij 75 = 22.0 ± 0.1 AB) for early-type galaxies at z = 0.83. This result shows that the 
completeness magnitude of our fields corresponds to ~ M* + 2 at the redshift of the pairs. 



3.1. Comparison with HHDFN and ACS-GOODS photometry 



Our approach to compute photometric redshifts (§4]) makes use of a training set with 
galaxies of known redshifts measured in the same photometric bands. Consequently, pho- 
tometric redshifts are very sensitive to small zero-point changes. In order to examine this 
point, we have made a comparison between our ph otometry with tho se available for the 
HHDFN (Hawaii Hubble Deep Field North) region JCapak et all l2004h and for the ACS- 
GOODS (Advance d Camera for Surveys - Great Observatories Origins Deep Survey) region 
rtCowie et al1l2004h . 



Although one region is contained in the other, the photometric and spectroscopic data 
available for them are different. A comparison with HHDFN is useful because it contains 
ACS-GOODS and has a photometric completeness similar to those of our fields. ACS- 
GOODS, on the other hand, has hundreds of mesured spectroscopic redshifts and will be 
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adopted as the training set for our photometric redshift method. Since the photometry of 
these fields (in UBVRIz') is diferent of our photometry (in g'r'i'z'), they were interpolated 
to the GMOS bands adopted here. 

First we have considered the HHDFN region which has photometry complete down to 
R = 24.5. We compared the photometric distribution of galaxies in our fields within the 
magnitude completeness limit of each band with the corresponding distributions using the 
interpolated magnitudes. We estimated the shift that is required in the zero-point of our 
photometric bands so that the median of the magnitude distribution (for galaxies brighther 
than the completeness limit) of HHDFN and ours match each other. We have used an 
iterative algorithm, adding to our magnitudes a shift obtained in each step until convergence. 
The median of the absolute value of the shifts is 0.07. 

Many are the possible sources of these shifts. A possibility is the magnitude interpolation 
required in this approach. Another one is cosmic variance, since the galaxy catalogues 
co nsidered here ar e not l arge enough. Indeed, a photometric study with SDSS data made 



by iFukugita et al.l (12004 ) shows that, besides Galaxy extinction, the principal cause for 
variations in number counts is the large-scale clustering of galaxies. This dispersion increases 
for smaller areas, being greater than 0.2 magnitudes for areas smaller than 0.01 deg 2 , as is 



our case. 



After appling the zero-point shifts in our magnitudes, we compared our galaxy number 
counts with those from the ACS-GOODS data in 0.5 mag intervals for objects brighter than 
the apperture magnitude z' = 22. It can be verified that, for all fields, we obtain a good 
match between our number counts and those of ACS-GOODS. 



4. Photometric redshift analysis 



Determining the redshifts of the galaxies in our fields allows us to separate the galaxies 
belonging to a possible cluster or group at the redshift of the pair from the foreground and/or 
background galaxies. Here we adopt photometric redshifts for this task. 

Photometric redshift estimation is often done by comparing the magnitudes of an object 
with the magnitud es of templates obtaine d with spectrop hotometry evolution m odels, as is 
the case of Zpeg (ILe Borgne et al.l |2002| ) and HyperZ (IBolzonella et al.l l2000l ). Here we 
adopt another approach: instead of a galaxy model, we use real data- magnitudes and 
spectroscopic redshifts- obtained in galaxy surveys. We compare the magnitudes of our 
galaxies with magnitudes in the same bands of real galaxies with known spectroscopic redshift 
to parametrize a local empirical relationship between magnitudes or colors with redshift. 
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This is done using a locally weighted regression algorithm (LWR) developed by our group 
(Santos et al. 2007, in preparation). The same type of data-driven approa ch is adopted in 
the ANNz photometric redshift package (IFirth. Lahav &: Somerville 1120031 ). which applies 
instead artificial neural networks for this task. 



4.1. Method 

LWR is an algorithm designed to provide a continuous non-linear mapping between sets 
of variables (e.g., Atkeson, Moore & Schall 1997). Our LWR method is discussed in detail 
and compared with other methods in Santos et al. (2007, in prep.). Here we only outline 
its main characteristics. The method works with magnitudes or colors (and even with other 
galaxy properties, like diameter or type) but here we use colors. 

The method works with two data sets: the training set (having known spectroscopic 
redshifts) and the test set (for which we want to calculate the redshift). Obviously, both 
sets must have the magnitudes and/or colors measured in the same bands and in the same 
photometric system. 

LWR establishes a linear relationship between colors and redshifts that is local because 
the redshift estimation in a given point in color space weights more heavily the data points 
in the neighborhood of this point than those more distant. The training set contains colors 
and spectroscopic redshift for all objects. From these values we build a redshift estimator 
which will be applied to our galaxies, in the test set. We assume that the local relation 
between colors and redshifts is linear: 

z(x) = a + a T .x = a + £™ =1 ajXj (2) 

where x is a vector containing the n colors of a given object, z(x) is the redshift and T stands 
for transpose matrix. For each object in the test set, with colors in point x, we determine the 
values of coefficients aQ...a n and then the redshift by minimizing the weighted \ 2 function 
with the N objects of the training set: 

X^Ef^Jfe-ao-a^) 2 (3) 

where u>j is the weight associated to the j—th data point. The locality of the fitting is assured 
by adopting a weight function which decreases as the euclidean distance d(x, x^-) between 
points x and Xj increases: 

/-c/ 2 (x,Xj)\ 
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The K parameter is a kernel-width that determines the "effective volume" around point x: 
only points within this sphere effectively affect the values of the parameters and the redshift 
estimate. This parameter was determined in this work by dividing randomly the training set 
objects in 2/3 for training (TS) and 1/3 for validation (VS). For galaxies in VS we computed 
Zphot with eq. [2] using galaxies in TS to obtain the coefficients. This procedure was then 
repeated one hundred times. For each realization of TS and VS, we compute the rms square 
deviation between z p hot and z spec for the objects in VS, a z , and choose as optimal K the 
value for which <j z is minimum. 

It is worth mentioning that redshift estimates with the LWR method are heavily depen- 
dent on the training set adopted (besides, of course, the set of colors available). In particular, 
the redshift accuracy increases with the size of the training set and depends strongly on the 
homogeneity of the photometric calibration of the training and test sets. 



4.2. Application to our sample 

We adopt in this work the (interpolated) photometric data and spectroscopic redshifts 
of the ACS-GOODS region as our training set. Since the method allows using colors or 
magnitudes for photometric redshift estimation, we have used colors (but none of the results 
reported in the next section depend of this choice). Using magnitudes we would have to 
limit our sample at z' = 22 (the ACS-GOODS spectroscopic completeness), but our sample 
photometry goes deeper; it is complete at least down to i' = 24. When we use colors, such 
limit is not necessary and we are able to estimate photometric redshifts for fainter objects, 
even without spectroscopic data for z' > 22. 

The value K = 0.33 was determined by the procedure described in the previous section 
as the median value of 100 simulations. The histogram in Figure [2] shows the redshift error 
distribution for all simulations. In what follows we consider the mean value, a z = 0.16, as the 
redshift error for the training set ACS-GOODS. Figure [2] also shows the comparison between 
z p hot and z spec for the 1/3 of galaxies from ACS-GOODS used for validation corresponding 
to the simulation with this mean value. 

Having obtained photometric redshifts for all fields containing quasar pairs, we may start 
looking for structures around the redshift of the pairs, that is, objects with z phot = z pair ±Az. 
We made experiments with Az equal to 0.1, 0.16 and 0.2, obtaining si milar result s . The refore 



we only present here results with Az = a z = 0.16. For comparison, iToft et al.l (120031 ). in a 
photometric study of the galaxy cluster MG2016+112 at z — 1, adopted Az = 0.25. Note 
that the error a z should depend of the photometric errors, the number of colors, and the size 
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of the training set. 



5. Results 

We present in Figure [3] the galaxy redshift distribution in the field of each quasar pair. 
All fields show a peak in the interval z G [z pair — <r z ,z pair + o z \. We now analyse some 
properties of the galaxy distribution in this redshift interval aiming at constraining the 
nature of the environment inhabited by the quasar pairs of our sample. 



5.1. Galaxy overdensities 

We must know the expected number of galaxies in this interval, to verify the significance 
of the galaxy excess around z pa i r . For this estimate we have assumed that the HHDFN 
region is representative of the overall galaxy distribution. This sample is appropriate for 
this analysis because its photometric depth is comparable to that of our fields. Photometric 
redshifts were obtained with the method discussed in the previous section. We then defined 
the galaxy overdensity in the interval z pa i r ± o z as 



n 



pair 



n H 



n H 



(5) 



where n pa i r and n# are the number densities of galaxies in this redshift interval for a given 
field and for the HHDFN, respectively. We have considered as galaxies in HHDFN all objects 
with z phot > 0. 

Values of 5 for each pair are shown in Table [3j Errors in 5 were determined assuming 
Poissonian errors for N pair and Nh- The overdensity 5 ranges from 0.6 to 1.6 and is significant 
in all cases. Note that these results are affected by cosmic variance, since we have used only 
one reference field and the area occupied by HHDFN (0.2 square degrees) is very small, 
so that Uh is affected by the galaxy clustering in the HHDFN region. We have arrived at 
similar results using the VIMOS VLT Deep Survey aroun d the Chandra Deep F ield South 
(ILe Fevre et al.l 12004 ) and the Gemini Deep Deep Survey ([Abraham et al. 1 120041 ) regions. 



5.2. Distribution of galaxies 



In order to investigate the clustering properties of the galaxies in the chosen redshift 
interval around a quasar pair, we calculated the median projected distance between galaxies 
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and compared them with the same quantities obtained with 1000 simulations of random 
uniform galaxy distributions with the same number of objects and the same projected area 
of the observed fields. We may then define a confidence level, CL, that a field presents a 
galaxy distribution more clustered than an uniform distribution: 

cl = (6) 

N s 

where N(A0 > Adf) is the number of simulated fields with median projected distances larger 
than that of the observed fields and N s is the total number of simulations. 

We summarize the results of this analysis in Table [3j Two pairs are strongly clus- 
tered (QP1355-0032, QP0110-0219), one is moderately clustered (QP1310+0007) and one 
(QP0114-3140) is not clustered at all. 



5.3. Richness 

In order to estimate the richness of o ur fields, we have adopted an approach similar 



to the traditional Abell's richness criterion (lAbell II1958T ). defined as the number of galaxies 



brighter than m 3 + 2 (where m 3 is the magnitude of the third brightest cluster member) 
within a radius of 1.5 Mpc of the cluster center. A cluster is considered rich if it 
contains more than 30 galaxies according to such a definition. 

The Abell's radius considering the cosmological model adopted here is 2.1 Mpc. Assum- 
ing that the brightest galaxy in the redshift interval z pair ± a z is the brightest cluster galaxy, 
we computed the number of galaxies brighter than i' 3 + 2 by scaling their number in each field 
to the Abell area (N esc = iV/E, where E = A par / ' AAbeii ~ 0.5). This result was corrected for 
contamination due to background/foreground objects using counts in the HHDFN region. 
The results are shown in Table [3j All but one of the putative clusters are rich, according 
with this criterion. The field of QP1310+0007 seems to be the poorest of our four fields, and 
is poor also with Abell's criterion. However, this is the pair with the brightest galaxy in the 
corresponding redshift interval among all quasar pairs, and its poorness may be an effect of 
galaxy counts, since they grow strongly for increasing magnitude. Furthermore, note that 
our fields are smaller than the Abell's radius, then the quasars could be in a poor cluster, 
group, or in the neighborhood of a cluster. On the other hand, the pair QP0114-3140, which 
is not rich by the results of Sections 5.1 and 5.2, has R=0 in Abell's classification. It is, 
then, appropriate to look for other richness estimators to confirm or not these results. 



Another useful richness indicator is the number of bright galaxies, assumed here as those 
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brighter than %'* + 1 present in the field. This number is estimated for galaxies in the pair 
redshift interval and is corrected with the corresponding HHDFN counts (scaled to the field 
area). The results are also presented in Table [31 All fields seems to contain a considerable 
number of bright galaxies. 



It is interesting to compare our results with those obtained by iPostman et al.l (120021 ). 
These authors studied a variety of Abell-like richness indicators in an I-band cluster survey. 
One of these indicators, Na^.5, is defined as the number of galaxies with magnitude between 
m 3 and m 3 + 2 within a radius of 666 h^l kpc. They show that this indicator is related to 
Abell's richness, Na, as Na,o.5 ~ OAANa- We have used this relation to estimate Abell's 
richness from Na,o.5- For 31 clusters with redshifts between 0.9 and 1.0, we obtain Na = 54 
galaxies. That means that fields have a richness R ~ 1, which may be compared with the 
numbers present in TableBl only QP1310 +0007 seems poorer than the clusters at comparable 
redshift studied by IPostman et al.l (120021 ) . 



5.4. The red sequence 

The red sequence is a characteristic of the color-magnitude diagrams of early-type galax- 
ies of groups and clusters. In a color- magnitude diagram these galaxies have very similar 
colors following a linear relation and their integrated colors are progressively bluer for weaker 
magnitudes. This relation is also known as color-magnitude relation (CMR). 

We have examined the red sequence in the (i f — z') x %' diagram of galaxies in the redshift 
interval of each pair. The use of the color {%' — z') is based on its capability to identify early- 
type galaxies, since at z ~ 1 the 4000 A break lies in the i' band and consequently the early- 
type galaxy color (i' — z') are very red. The fields around quasar pairs QP1310+0007 and 
QP01 10-0219 (Figures H] and [6] - top-right) present a peak in the color distribution at 0.6 < 
i'—z' < 1.0. Comparing our data with a similar distribution for HHDFN, we verify that these 
peaks represent an excess of 1.7 a and 3.3 a, respectively. Therefore, in this interval, we would 
expect to find a red sequence in the color-magnitude diagram. Indeed, for QP0110-0219 we 
note clearly that the galaxies fo r m a r ed sequence in %' — z' ~ 0.8 (Figure [6]- top-left), the 



value obtained by iTanaka et al.l (120061 ) in spectroscopically confirmed structures at z ~ 0.9. 
Besides, if we consider the projected distribution of these red-sequence galaxies (Figure [6]- 
bottom), we notice that they have a filamentary-like distribution similar to what is observed 
in other z ~ 1 clusters, and considered typical of clusters in process of formation (e.g., Toft 
et al. 2003). 



The red galaxies of QP1310+0007 present a broad distribution in the color-magnitude 
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diagram (Figure H]- top-left). They also present a clump-like projected distribution (Figure 
H]- bottom). The other two fields have less-significant red sequences (Figures [5] and [71). The 



cluste r CL1604+4321, at z ~ 0.9, the less massive of the clusters studied by iHomeier et al. 



(120061 ). presents a lack of bright elliptical galaxies (~ M*). The authors suggest that this 
cluster has not yet had time to complete the red sequence. This may be also the case for 
the structures associated to QP1355-0032 and QP0114-3140. 



5.5. Properties of the fields around quasar pairs 

The properties of the environment associated with each quasar pair are summarized in 
Table HI We now discuss each pair individually. 



5.5.1. QP1310+0007 

This quasar pair is formed by J131046+0006 (a radio-quiet object) and J131055+0008 (a 
radio-loud quasar) at redshifts 0.925 and 0.933, respectively. They have an angular separa- 
tion of 177 arcsec, corresponding to 1.4 Mpc in the adopted cosmology. Its density contrast 
is the smallest among all quasar pairs. However, the galaxy distribution analysis shows that 
the galaxies in this field are clustered at some degree, i.e., the median projected distance 
between galaxies is smaller than in a random uniform field in 67 % of the simulations. This 
field has been classified as poor with the Abell's criterion, and we have found 20 galax- 
ies with magnitude i' < i'* + 1. The galaxy color distribution shows a prominent peak in 
0.6 < i' — z' < 1.0, corresponding to the red sequence. The red galaxies present a clumpy 
distribution, but without central condensation. The presence of a significant amount of 
early-type galaxies plus the relative poorness of this field indicate that this can be the seed 
of a structure that can became a rich galaxy cluster at z — 0. 



5.5.2. QF '1355-0032 



J135457-0034 and J135504-0030 constitute this pair; the first is radio-loud and the 
second is radio-quiet. The projected separation between them is 252 arcsec, or 2.0 Mpc, 
with redshifts 0.932 and 0.934, respectively. The redshift interval z pair ±a z shows the largest 
galaxy excess among all quasar pairs discussed in this work. The median projected distance 
between galaxies resulted smaller than in a random uniform field in 98.5 % of the cases, 
meaning that these galaxias are strongly clustered. This is the richest field in our sample 
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accordingly to Abell's criterion and also the one with the largest number of bright galaxies 
(i' < %'* + 1), however its red sequence is modest and the red galaxies do not present a 
clustered distribution. Its richness and number of bright galaxies are the major indications 
that this quasar pair is probably in a galaxy cluster. 

5.5.3. QPO '11 0-0219 

This quasar pair is formed by a radio-loud quasar (Q 0107-0235) and a radio-quiet 
quasar (PB 6291) at redshifts 0.958 and 0.956, respectively. The angular separation of 77 
arcsec (0.6 Mpc) is the smallest of the sample. The overdensity in the redshift interval is 
significant as for the other pairs. The galaxy distribution is the most clustered of all samples, 
accordingly to the CL values in Table [3j The field is rich by Abell's criterion, but presents 
only 12 bright galaxies. The red sequence is clearly present in the color-magnitude diagram 
at %' — z' ~ 0.8. The red galaxies present a filamentary-like distribution and there is a galaxy 
excess around the radio- loud quasar. These results indicate that QP01 10-0219 is indeed a 
rich cluster. 

Moreover, we have verified that QPO 11 0-02 19 has been serendipitously detected (but 
unreported) in X-ray with a pointed ROSAT PSPC observation of 6.6 ks. We have estimated 
the bolometric X-ray luminosity assuming that all detected flux (background corrected) 
comes from the ICM: Lx,bo\ ~ 5 x 10 45 ergs s _1 . Such luminosity is well above a typical 
cluster X-ray luminosity and may be contaminated by the X-ray emission from one or both 
quasars. On the other hand, the typical quasar X-ray luminosity [2-10 keV] is around 
10 44 erg s _1 , thus the quasars in the pair may not account for all X-ray emission. Besides, 
the total emission [0.5-8.0 keV] within 3 arcmin is about 10 times higher than the typical 
quasar emission in this band; therefore the observed X-ray flux is consistent with emission 
from quasars and a possible cluster around them. This is the single pair in our sample 
detected in X-rays so far. The other fields were not detected in the Rosat All Sky Survey, 
nor in the pointed observations. 

5.5.4. QP0114-3140 

This pair is formed by radio-quiet quasars. J011441-3139 has z = 0.974 and J011446- 
3141 has z = 0.968. The separation between them is 144 arcsec (1.1 Mpc). This field shows 
a significant overdensity, no clustering, and no red sequence. It also seems rich with Abell's 
criterion and has a comparatively large number of bright galaxies. The evidence for the 
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presence of a rich cluster at the redshift of the quasar pair is not as compelling in this case, 
compared with the other 3 pairs. 

6. Summary 

We have studied the environment traced by quasar pairs at z ~ 1, using images in 
g', r', i', and z' bands obtained with GMOS at Gemini North and South. In order to 
identify galaxies in a redshift interval close to that of the quasar pairs, we have estimated 
photometric redshifts with the LWR method, using ACS-GOODS data as a training set. 
The rms dispersion of the difference between our photometric redshift and the spectroscopic 
redshift in the training set is a z = 0.16. We have adopted the interval z = z pair ± a z for the 
analysis of the pair environment. 

When compared with the HHDFN region, all fields show a significant overdensity in the 
redshift interval of the pair. In all cases this excess is larger than 3.5 a. 

We investigated the clustering of the galaxies near the pair by estimating a confidence 
level, CL, that the galaxies are more concentrated than in a uniform distribution. We have 
also estimated the richness of each redshift interval with a variant of Abell's criterion, as 
well as by the number of bright galaxies. We verified whether a red sequence is present and 
the form of the projected distribution of red galaxies. 

The analysis indicates that probably three out of our four quasar pairs are members 
of galaxy custers. For one of the pairs we did not find strong evidence for it: QP0114- 
3140 could be in a poor cluster, group, or in the neighborhood of a cluster, since our fields 
are lesser than Abell's radius. Taken at face value, this result shows that quasar pairs are 
indeed good tracers of the large scale structure at high z. However, with only four quasar 
pairs in our sample we are not able to say at what level targeting a quasar pair increases the 
probability of finding a rich galaxy cluster as compared to targeting a single quasar. A study 
of larger and homogeneous samples would be necessary to clarify this point. An extension 
of our work to other redshifts may be also useful and may provide interesting clues on the 
evolution of large-scale structure and galaxy clustering. 

This work is based on observations obtained at the Gemini Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partnership: the National Science Founda- 
tion (United States), the Particle Physics and Astronomy Research Council (United King- 
dom), the National Research Council (Canada), CONICYT (Chile), the Australian Research 
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Council (Australia), CNPq (Brazil) and CONICET (Argentina). 
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Fig. 1. — Star/Galaxy separation for the pair QP0110-0219. Stars represent objects with 
FWHM < seeing and filled circles represent objects with FWHM > seeing. We have consid- 
ered as galaxies all objects with class-star < 0.9 (dashed line). 
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Fig. 2. — Left: photometric redshift error for 100 simulations. Right: comparison between 
photometric and spectroscopic redshifts for the 1/3 of galaxies from ACS-GOODS used 
for validation, corresponding to the simulation with photometric redshift error equal to 
<t 2 = 0.16; the continuous line is the equality line and the dashed lines correspond to 1 and 
3 a z . 
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Fig. 3. — Distribution of photometric redshifts for each field containing a quasar pair: a) 
QP1310+0007, b) QP1355-0032, c) QP0110-0219, and d) QP0114-3140. The hatched area 
in each panel corresponds to z pair ± a z . 
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Fig. 4. — Results for the pair QP1310+0007. Galaxies within z pair ± a z . Top-left: color- 
magnitude diagram; galaxies in the red sequence are represented as filled circles. Top-right: 
histogram of the color distribution. Bottom: projected distribution of galaxies. The quasar 
pair is represented by stars and the galaxies in the red sequence as filled circles. 
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Fig. 5.— Same as Fig. Shut for QP1355-0032. 
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Fig. 6.— Same as Fig. 01 but for QP0110-0219. 
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Fig. 7.— Same as Fig. Sbut for QP0114-3140. 
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Table 1. Sample Characteristics 



Quasar Names 




a 




<5 




z 


A6» 


Quasar Pair Name 




(■ 


12000) 


(J2000) 




(arcsec) 




J131046+0006* 


13 


10 


46.2 


00 06 


33 


0.925 


177 


QP1310+0007 


J131055+0008 


13 


10 


55.9 


00 08 


14 


0.933 






J135457-0034 


13 


54 


57.2 


-00 34 


06 


0.932 


252 


QP1355-0032 


J135504-0030* 


13 


55 


04.7 


-00 30 


20 


0.934 






Q 0107-0235 


01 


10 


13.2 


-02 19 


53 


0.958 


77 


QP0110-0219 


PB 6291* 


01 


10 


16.3 


-02 18 


51 


0.956 






J011441-3139* 


01 


14 


41.8 


-31 39 


25 


0.974 


144 


QP0114-3140 


J011446-3141* 


01 


14 


46.4 


-31 41 


31 


0.968 







* Radio-quiet quasars 
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Table 2. Observations 



Pair 


Telescope 


<?' 


texp (seconds) 
r' i' 


z' 


Identification Number 


QP1310+0007 


Gemini N 


9 x 300.0 


6 x 200.0 


11 x 350.0 


8 x 450.0 


GN-2003A-Q-2/GN-2005A-Q-19 


QP1355-0032 


Gemini N 


13 x 300.0 


6 x 200.0 


6 x 350.0 


7 x 450.0 


GN-2003A-Q-2/GN-2005A-Q-19 


QP0110-0219 


Gemini N 


10 x 300.0 


6 x 200.0 


8 x 350.0 


8 x 410.0 


GN-2003B-Q-2/GN-2004B-Q-24 


QP0114-3140 


Gemini S 


7 x 300.5 


6 x 200.5 


7 x 350.5 


7 x 410.5 


GS-2003B-Q-3/GS-2004B-Q-17 
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Table 3. Clustering Properties 



Pair 


8 


median 

(arcmin) 


C 'median 

(%) 




N(i' <i' 3 + 2) 


N esc (i' < i' 3 + 2) 


AT(i' < i'* + 1) 


QP1310+0007 


0.58 ± 0.14 


2.7 


67.0 


20.35 


6 (R < 0) 


13 (R < 0) 


20 


QP1355-0032 


1.59 ± 0.19 


2.6 


98.5 


21.06 


95 (R = 2) 


203 (R = 4) 


58 


QP0110-0219 


0.70 ± 0.14 


2.4 


100.0 


21.29 


35 (R = 0) 


72 (R = 1) 


12 


QP0114-3140 


0.86 ± 0.23 


2.8 


0.5 


20.63 


34 (R = 0) 


95 (R = 2) 


36 
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Table 4. Summary of the quasar pair properties 



Pair 


5 


CL 


N A 


RCM 


C/F a 


X-rays 


QP1310+0007 


ok 


ok 


X 


ok 


ok 




QP1355-0032 


ok 


ok 


ok 


X 


X 




QP0110-0219 


ok 


ok 


ok 


ok 


ok 


ok 


QP0114-3140 


ok 


X 


ok 


X 


X 





a Clustcr-likc or filament-like distribution. 



